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УДК 631.362.001.573 
Ю.И. ЕРМОЛЬЕВ 
МОДЕЛИРОВАНИЕ ВОЗДУШНО-РЕШЕТНОЙ ЗЕРНООЧИСТИТЕЛЬНОЙ МАШИНЫ 


Приведена математическая модель процесса функционирования воздушно-решетной зерноочистительной 
машины, показаны условия и результаты этого моделирования. 
Ключевые слова: процесс, показатели сепарации, математическая модель, результаты моделирования. 


Введение. Одним из путей первичного проектирования воздушно-решетных зерноочистительных 
машин можно считать их параметрический и структурный синтез. Приведен пример параметриче- 
ского синтеза машины с использованием системного моделирования. 

Построение математической модели. Математическую модель процесса функционирования 
воздушно-решетной зерноочистительной машины (ВРЗОМ) как замкнутой квазистатичной системы 
с заданной Аи -й функциональной схемой в общем виде можно записать [1]: 


Ефём = |8, А, С, [м (х),Тзм ()|-> тах, (1) 
Ам а Пей, хЕСи (х,и); (2) 
Ь < |, ра бы 9 9 (3) 
Здесь Р — вектор входных воздействий на принятую в зерноочистительной машине систему опе- 
рации. 
Е=,а,„И,М(, } 5, } Уо(в), о (Н) }, (4) 
А - вектор управляющих факторов системы, обеспечивающих ее функционирование, 
д = Вий, и» Г, (Ви), .КВп),Р. (В), = (Н), Хо» (В).Кам (®} (5) 


Векторы входных и управляющих воздействий на пневмосепаратор 1 ГК и Ад, на решет- 


ный модуль — Аь и Аь, на пневмосепаратор 2 — Ел. и Ан, показаны на рис.1. 
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Рис.1. Структурная схема воздушно-решетной зерноочистительной машины 
с многоярусными решетными модулями 
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В формуле (1) х - набор технологических операций из С’, (х,и) множества, реализуемых 
в К,, функциональной схеме ВРЗОМ. С, К (х), 7, (х) — функционал, определяющий показатели 


технологического процесса в ВРЗОМ (полнота выделения }-х компонентов, потери зерна, содер- 
жание в очищенном зерне других }х компонентов), для принятой К’, (х) схемы. 


Выходные показатели функционирования зерноочистительной машины определяются век- 
тором Вм, независимые аргументы которого случайные в вероятностно-статистическом смысле 
величины: 


йе 


Вы — Евы, ВВ, бы, бз ОО: 
Пе оба } 
где А, -— критерий эффективности реализации технологического процесса ВРЗОМ; 5 „- Полнота 


(6) 


выделения из зернового материала отделяемых /-х компонентов; [ О. — потери зерна на 
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хр операции и общие; О,,О,,О„ - масса выхода очищенного зерна, фуражных отходов и 
отходов разных категорий; Упл› Ух - Масса }-х и всех компонентов, содержащихся в зер- 
новом материале после реализации х/-й операции; у Ух — Масса у-х и всех компонен- 


тов, выделенных из зернового материала после реализации хт-йЙ операции; Р,.. - содержа- 


пря 
ние 7-ГО компонента в массе очищаемого материала и в выделяемых фракциях Бут; 2 Ее 
полнота прохода Ко) компонента в очищаемый зерновой материал после выполнения хх -Й 
операции. 


Еф =1- ^^, (7) 
>.а, 
2 
функция цели: 
Е —> тах. (8) 


Обозначения элементов, входящих в выражения (1)-(8), определены в предыдущей работе 
автора [1]. 

В общем виде плотность вероятности распределения зернового материала по ширине В 
пневмоканала и решетного модуля — Ло (В), а плотность вероятности распределения скорости 


воздушного потока в пневмоканале — Г, (В,). Тогда количество зернового материала О, попа- 
дающего на каждый р-й участок (р=1, 2, ..., е) К-й ширины рабочего органа, 


, | о(вдав 
ре (9) 
| о(в)ав 


а средняя скорость воздушного потока на р-м участке. 


| х.(вуав 
Г. = т | (10) 
| х.(вуав 
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Полнота выхода 7-х компонентов зернового материала в очищенное пневмосепаратором 





зерно: 
е 
хо в’ (О Г 
| ВЕ. (11) 
пА7 
У Ое 
где 4/(О Г „) — полнота выхода 7-ГО компонента в очищенное пневмосепаратором зерно для ус- 


ЛОВИЙ сепарации, определенных подачей в зернового материала и скоростью , воздушно- 
го потока на р-м участке пневмосепаратора 
(о — )= Ло а ›М ,Л.(В).У (В.В ›й на (12) 
У) р Р р 70 У у О шо ш т 
где М»; — математические ожидания средних скоростей витания 7-Х компонентов в исходном мате- 
риале, Ва; — ширина 1-ГО пневмоканала, Г — глубина -ГО пневмоканала, :-— условия ввода 
зернового материала в пневмоканал. 


Конкретные математические модели, определяющие & для вертикальных пневмокана- 


п) 
лов, пневмосепараторов с активным вводом зернового материала в вертикальный или наклонный 
пневмоканал определены [1]. 

При выполнении ограничений (1)-(3) можно рассчитать показатели процесса сепарации: 

— содержание сорных примесей (;=1, 2, ..., с) в очищенном зерне 


66 =>, а,=, | >. авы (13) 


— полнота выделения &». пневмосепаратором в отходы сорных (1=1, 2, ..., с) примесей 


бьс - У а-еы а, р (14) 


Ее 
— содержание гх компонентов в очищенном зерне 


Ь 
Ь, =а,Ес;/ уз Я пк; |, (15) 
1 
— полнота выделения то компонента зернового материала в пневмосепараторе 
&, =(-&); (16) 
— содержание то компонента зернового материала в легких отходах 
Ь 
6, =а,(1-в)/| У а, 4-2.) (17) 
у 


— масса отходов и масса очищенного в пневмосепараторе зерна 
Ь 


О, =» Оа,1-=шу), КГ/(м-с); Ок =О- Оо, кг/(м.<). (18) 
= 
Построенная математическая модель процесса сепарации гетерогенной сыпучей среды 
адекватна [2] и позволяет, используя известные методы параметрического и структурного синте- 
за, проводить многомерный анализ процесса и оценивать рациональные Параметры пневмосепа- 
ратора для задаваемых условий и ограничений на показатели его функционирования. 


Входное воздействие на рассматриваемый решетный модуль выражается вектором Р,, 


независимые аргументы которого случайные в вероятностно-статистическом смысле величины. 
г” И 
19 =. | „аз, М, } в (ь.) (в), (Н.Г) ь (19) 
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Активные средства, определяющие эффективность функционирования решетного моду- 
ля, найдем вектором А, 


А. =@,,В,,А,.м,, т... Вы,с, К, (х)|, (20) 


где все показатели определены [1]. 

Различные базовые решетные модули, которые можно использовать в ВРЗОМ, математи- 
ческие модели, описывающие процессы функционирования базовых и комплексных (составлен- 
ных из различных структур базовых модулей) решетных модулей, определены [1]. 

Рассмотрим, для примера, в общем случае многоярусный решетный модуль (рис.2) с неод- 
нородными по длине (Ду,.../„) и высоте (1, 2, ..., К) решетами и с различными кинематическими 


параметрами. В общем виде на ярусы решет поступает 41 ;,4> ;,....4„;›4х; Количество /-го компо- 


нента сыпучего материала с их торцов, определяемое плотностью вероятности распределения 
Ло(Н ) случайной величины О’ (подача сыпучего материала на ярусы (1, 2, ..., К решет)) по вы- 


соте Н расстановки ярусов решет: О=О, (18), содержание -х компонентов в сыпучем материале 
Ь; определяется из выражения (15). Распределение случайной величины а; по ширине Вр се- 
параторов в ярусах определяется плотностью вероятности /„, (В) В качестве допущения здесь и 


далее будем считать поступление сыпучих материалов на ярусы решет постоянными по времени 
(процесс квазистатичный) и с равномерным распределением у-х компонентов в сыпучем мате- 


риале (7 =1, 2, ...). 








и 
Чар “# 924 та а к ты 


я 


ук "Я а ны о т 85; 
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Улк}  'Упк? = ыы 
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Рис.2. Многоярусный решетный модуль (а). Плотность Ло (НН) вероятности поступления зернового материала 


на А-е решетные ярусы и плотность вероятности Леа (В') его распределения по их ширине В, (6) 


Полноту просеивания 7-го компонента на 6-м решетном ярусе можно определить по из- 
вестному выражению [1]: 


1 
> Чъа,е5 (,) 


=5: И (21) 
Ча 1 
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де 4» - приведенная подача сыпучего материала на р-й участок (по ширине В) 


В 
(р=1,2,....е) п-го решетного яруса, Ч бр = Во: к=В|е — ширина р -го участка яруса 


решет; Ко, - вероятность попадания случайной величины О на р-й участок по ширине 
д -го решетного яруса. 


Бор = Ро Рбра . (22) 


Ор 
Здесь Ро - вероятность попадания случайной величины О на ге: (2) участок (6-й 


решетный ярус): 


Гони [вв 


В. В 


ра — 


О (23) 


о [убнин Глав. 


При этом подача 4, сыпучего материала на б-й решетный ярус 
45 = ОРлок; (24) 
55; (4. ) — полнота просеивания /-го компонента на р-м участке 6-го решетного яруса при по- 
даче 4, сыпучего материала на этот участок, определяется из различных выражений в зависи- 
мости от многих факторов: 
55 (Ч |= Л(О,а „у, ина Ви» Ви, Миьп), 
где /,,у — вид сыпучего материала (зерно, семена трав и др.) и его плотность; Г, — длина т- 
го решета в -го решетного яруса; 4;, — рабочий размер отверстий т -го решета д -го яру- 
са; „ВАМ, - Угол наклона к горизонту, направленность, амплитуда и частота ко- 
лебаний т -го решета в д-м ярусе; и -— количество решет в ярусе. 


[ Го(Н)ан | (В)авВ 





ро Г 45 5р (Ч р) 
7" | (ван | и (вав 
вы = и=1 р=1 | (25) 
| Ч5 а) 


И 
ГДЕ 55 = („пор а 
При известной полноте просеивания =, /-го компонента на отдельном ярусе решет 


(9=1,2...., К) полнота просеивания /-го компонента на всех ярусах решетного модуля определя- 
ется из выражения 


|. 
> Ча 55} 
__ 0=1 
| о иииетвИ (26) 
Оа;к 
где а’ - содержание /-го компонента в сыпучем материале, поступающем на 6-й решетный 


ярус. 
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Сход /-го компонента с р -го участка д„-го решета определяется из выражения 


и — 25, В. р 
К =аза,] | а, „(® }ехр ОЕ ЕЕ, ти 8 (27) 
1=1 бр 


Принимая, что для каждого р -го участка и для каждого бт -го решета параметры сепа- 
рации аз»; ( ), Би, Ги; Различны, но постоянны в пределах ограниченной длины бт -го решета, 


массу просеивания /-го компонента сыпучего материала на длине Х р -го участка бт -го ре- 
шета можно определить из выражения 


7 275" Вирт Изт-Х 
7 бт) (Хх) = 15 (т-1)} Е Ч 5 ту (7 )ехр г —_—_—_—_ / (28) 


ь уз 6(т-1) 
7 


а полноту просеивания /-го компонента на длине / бт -го решета модуля - из выражения 


/ 
>. Ч бтр Я бт б6т] (Ч бт > Г) 





ООВ (29) 
[оапан | 
О Па-=5 а! 
) | осан "^ 
6=1 


Полнота просеивания /-го компонента =„, через все т -е решета А ярусов определяется 


из выражения 


| хадан 
вт = Фовы —- 


м [ан 


6=1 


(30) 


Полнота просеивания /-го компонента =, через все решета ^ -ярусного решетного моду- 

ЛЯ 
8; =2(и1)о; т | > Не би; (31) 
где 8(и-1)0; — полнота просевания 7-го компонента через все (п-1) последовательные решета в 


К ярусах. 
Полнота схода &сх»„, Л-го компонента со всего от -го решета 


Г 
>. Ч5"рЧ ;8 сХот] @ы, 
р= 


Ч5т@ } 
а со всех т-х решет А -ярусного решетного модуля 
А 
бсхт = _ _ (33) 
5=1 
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Полученные выражения (19)-(33) позволяют оценить проход /-х 4„з„; И всех 4,„„ компо- 


нентов сыпучего материала через би -е решета отдельных ярусов и всех А ярусов О 


пот] ! пот 


и 
содержание /-х компонентов в проходе каждого т -го решета 6 -го яруса Б;„, и всех т -х ярусов 
В 

Входное воздействие на второй пневмоканал определяется вектором т ‚ активные сред- 


ства — вектором Я (рис.1). 


Математическая модель оценки показателей функционирования второго пневмоканала 
аналогична (9)-(12) при подаче зернового материала в пневмосепаратор 2, равной выходу очи- 
щенного в решетном модуле зерна. Проведенные исследования показали, что построенные мате- 
матические модели с 95-процентной (полнота выделения сорных примесей — 90%) доверительной 
вероятностью адекватно описывают процесс сепарации в воздушно-решетной зерноочиститель- 
ной машине [1]. 

Параметрический синтез. Параметрический синтез функциональной схемы и задаваемой тех- 


ническим заданием вариации величин аргументов векторов входных Е (4) и управляющих А (5) 
воздействий проводится с использованием соответствующих математических моделей методом 
нелинейного программирования — методом сканирования с ограничениями, так как отсутствует 
унимодальность функции цели (8) и нелинейные ограничения (2)-(3). 

На втором иерархическом уровне синтезируется функциональная подсистема решетного 
модуля (РМ) с делителем зернового потока (ДЗП) для ВРЗОМ. 

Для синтезированных рациональных параметров рассматриваемой структуры РМ (напри- 
мер, схема РМ, число решетных полотен, форма и размеры отверстий решет, кинематические па- 
раметры: частота, амплитуда колебаний, угол наклона решет к горизонту, угол направленности 
их колебаний) проводится многомерный анализ процесса функционирования РМ при различных 


величинах вектора выходных РГ, (19) и управляющих А, (20) воздействий, определяемых выход- 


ными показателями функционирования ПС. 

Эффективность сепарации зерновых (пшеница, рис) и семян подсолнечника с учетом ог- 
раничений размеров решетных полотен (ТУ 23.2.2068-89), изготавливаемых в России, а также 
приводов решетных модулей и систем их очистки от забивания отверстий выявлена для одно- 
ярусных трехрешетных модулей. 

Определен рациональный решетный модуль с параллельной и последовательной схемами 
сепарации зернового материала, функциональные показатели которых соответствуют заданным, 
для очистки зерна с учетом проведенного параметрического и структурного синтеза комплексных 
решетных модулей [1] 

С использованием математической модели (1)-(33) был разработан алгоритм программ 
для ЭВМ ПК, позволяющий проводить параметрический синтез (8) с учетом принятых ограниче- 


ний (2)-(3) и значений аргументов векторов ГР, 4, А,,В„ [1]. 
Для проведения многомерного анализа оценки основных показателей функционирования 


обоснованного варианта ВРЗОМ (рис.3) при вариации аргументов вектора РК (рис.1) осуществле- 
но моделирование процесса системной сепарации зерна. 
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Условные обозначения: 


——3» - ввод обрабатываемого материала; 
— основной поток зернового материала; 








— В — — воздушный поток; 

_ .— - подсев; 

— „— - легкие примеси с воздухом; 
-$— -Фураж; 

—/— — легкие примеси; 

— х—  — крупные примеси; 


| 


В] 
$ $ ++ 


4 х— 
ч—- 


—»№ выходы фракций 


Рис.3. Технологическая схема воздушно-решетной зерноочистительной машины: 
1- камера приемная; 2- скельператор барабанный; 3, 4-— каналы аспирационные; 
5 — камера осадочная; 6 - 4-ярусный решетный модуль 


При моделировании приняты следующие допущения: плотность распределения вероятно- 
сти поступления зернового материала по ширине решетного модуля подчиняется нормальному 
закону с коэффициентом вариации 20,26%; плотность вероятности распределения зернового ма- 
териала по высоте решетных ярусов характеризуется коэффициентом вариации г. ), изме- 


няющимся в пределах от 0 до 54,4%. 

Результаты моделирования. Выявлено (рис.4,а), что в рассмотренном интервале подач изме- 
нение коэффициента вариации неравномерности поступления зернового материала на отдельные 
решетные ярусы несущественно влияет на выделение из него крупных сорных примесей и дроб- 
леной соломы. 


Полнота выделения мелких сорных примесей =, (рис.4,6) в интервале подач от 2,775 до 


4,635 кг/с незначительно уменьшается при изменении коэффициента вариации от 0 до 54,41%. 
Для интервала подач 4,635-8,58 кг/с =, незначительно снижается при изменении коэффициента 


вариации от 0 до 23%. При дальнейшем увеличении ло(Н ) происходит существенное снижение 


полноты выделения мелких сорных зерновых примесей. Агропоказатель Б.<1% (Б. - содержа- 


ние сорных примесей в очищенном зерне) выполняется на всем рассматриваемом диапазоне ва- 
риации подач О и Хо(Н) (рис.5, а). 


Существенное влияние изменения коэффициента Хо(Н ) и подачи оказывают на выделе- 


ние из рассматриваемого материала зерновых примесей = „. на всем интервале изменения подач. 


З 


Например, для О =3,705 кг/с при изменении Хо(Н ) от 0 до 54,4% =„„ снижается на 5%, а для 


З 


0 =8,115 кг/с при изменении АН ) в том же интервале =„. снижается на 7%. При этом агро- 


З 


показатель Р.< 2% (содержание зерновых примесей в очищенном зерне) выполняется при раз- 


личных значениях Хо(Н ) для разных подач. Так, при ло(Н )=24% максимально допустимая по- 
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дача составляет 6,495 кг/с (23,2 т/ч). Для О =5,65 кг/с и более увеличение Хо(Н ) существенно 
влияет на потери зерна в фуражные и неиспользуемые отходы. Так, для О = 6,495 кг/с макси- 


мально допустимое значение Хо(Н ) составляет 35,66%, а для О=8,115 кг/с Хо(Н ) „=181%. 


тах 





17,66 Л, % 


а) 6) 


Рис.4. Зависимость полноты выделения компонентов из зернового материала, очищенного воздушно-решетной машиной 
от подачи О’ зернового материала в машину и величины коэффициента вариации Л, неравномерности поступления зер- 


нового материала на отдельные решетные ярусы: а- крупных „кн и 6- мелких сорных 5, 
см 


Выявлена усредненная допустимая неравномерность поступления зернового материала на 
решетные ярусы в пределах Хо(Н )=28-30%. При синтезе основных параметров зерноочисти- 


тельной машины для различных основных величин О и / (Н ) усредненные рациональные па- 


раметры: Г, =6 м/с, и = 7,92 об/с, В, =1,7 мм, Ь› =2,0 мм, ВБ, = 3,6 мм. 


Ш 07-075 
9055-97 
008-055 
Ш055-05 
005-055 


Я по’ % 





а) 6) 


Рис.5. Зависимость содержания в зерне, очищенном воздушно-решетной машиной, сорных примесей Б, (а) и чистоты 


Я (6) зерна от подачи () зернового материала в машину и величины коэффициента вариации Л, неравномерности 


поступления зернового материала на отдельные решетные ярусы 
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Технические науки 


Выводы. Проведенное моделирование с использованием методов параметрического синтеза по- 
зволило проанализировать количественные зависимости показателей функционирования воздуш- 
но-решетной зерноочистительной машины при широкой вариации факторов, воздействующих на 
этот процесс для принятых условий и технологических ограничений, и синтезировать основные 
параметры ее работы [1]. 
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